Inter-isotope determination of ultracold rubidium interactions from three 

high-precision experiments 
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Combining the measured binding energies of four of the most weakly bound rovibrational levels 
of the 87 Rb2 molecule with the results of two other recent high-precision rubidium experiments, 
we obtain exceptionally strong constraints on the atomic interaction parameters in a highly model 
independent analysis. The comparison of 85 Rb and 87 Rb data, where the two isotopes are related by 
a mass scaling procedure, plays a crucial role. Using the consistent picture of the interactions that 
thus arises we are led to predictions for scattering lengths, clock shifts, Feshbach resonance fields 
and widths with an unprecedented level of accuracy. To demonstrate this, we predict two Feshbach 
resonances in mixed-spin scattering channels at easily accessible magnetic field strengths, which we 
expect to play a role in the damping of coherent spin oscillations. 
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After the first realization of Bose-Einstein condensa- 
tion (BEC) in a dilute ultracold gas of rubidium atoms 
|jj , experiments with the two isotopes 87 Rb and 85 Rb fur- 
ther lead to an amazingly rich variety of BEC phenom- 
ena, ranging from the controlled collapse of a condensate 
with tunable attractive interactions B to the realization 
of an atomic matter wave on a microchip || . Because of 
the large number of groups that have started doing exper- 
iments with these atomic species and the growing com- 
plexity and subtlety of the planned experiments, there 
is a clear need for a more precise knowledge of the in- 
teractions between ultracold rubidium atoms in the elec- 
tronic ground state, since these determine most of the 
properties of the condensate. For instance, despite a 
widespread interest, until now to our knowledge no ex- 
perimental group has been able to locate the predicted 
Pj magnetic-field induced Feshbach resonances that can 
be used to tune the interactions between ultracold 87 Rb 
atoms. Being able to switch on or off these interactions 
at will by a mere change of magnetic field may well be 
one of the main assets of matter waves compared to light 
waves in the new matter wave devices. In an atomic 
interferometry device, in particular, a nonlinear inter- 
action between interfering waves may be introduced or 
eliminated by changing a field applied at the intersection 
point. 

In this Letter, combining the results of three very re- 
cent high-precision observations, we come close to a com- 
plete and model-independent specification of the inter- 
action properties of ultracold rubidium atoms. The fact 
that two isotopes 85 Rb and 87 Rb are involved in the mea- 
surements makes the constraints exceptionally strong and 
also increases the predictive power: the interaction prop- 
erties of any other fermionic or bosonic isotope with mass 
number 82, 83, 84, or 86 are now known with about the 
same precision. Using mass scaling to relate the different 
isotopes we are able for the first time to deduce for each 



of the isotopes the exact numbers of bound Rb2 states 
with total spin S = (singlet) and 1 (triplet). As an 
illustration of the predictive power we predict 

two Feshbach resonances in mixed-spin scattering 
channels for 87 Rb at easily accessible fields that could 
lead to new time dependent phenomena in coherent spin 
oscillations and spin waves. There are numerous effects, 
such as spinor condensate energy differences, which are 
proportional to differences of scattering lengths. Because 
these differences are unusually small in Rb, the potentials 
must be very accurate to calculate them to reasonable ac- 
curacy. 

The first of the three high-precision experiments is the 
recent measurement of four of the highest bound rovibra- 
tional levels of the 87 Rb2 molecule with 10 kHz precision 
j5|. The second experiment is the improved characteri- 
zation H of the elastic scattering near a Feshbach reso- 
nance in 85 Rb, leading to a more precise determination of 
the resonance field Bq = 154.9(4) G and the nearby field 
strength B' = B + A = 165.85(5) G, where the scatter- 
ing length goes through zero (A is the (elastic) resonance 
width). The third experimental ingredient going into our 
analysis is the measurement JtJ of 12148 transition fre- 
quencies between X X E+ vibrational levels of the ( 85 Rb)2, 
( 87 Rb)2, and 85 Rb 87 Rb molecules, leading to a highly ac- 
curate singlet Rb + Rb potential ||. Moreover, within 
the accuracy of this experiment a comparison of levels 
for the three studied isotopomers shows no sign of Born- 
Oppenheimer break-down effects, i.e., the observed levels 
agree with a simple radial Schrodinger equation contain- 
ing a common singlet potential V<j(r) and the reduced 
atomic mass. 

This set of extremely precise measurements calls for 
a very careful construction of the interatomic total spin 
S = and 1 potentials, depending on the interatomic 
separation r. We combine the singlet potential of Ref. Q| 
with a long-range part equal to the difference Vdi sp —V exc h 
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of a dispersion term and an exchange term, starting at 
a variable radius r$ between 21 and 23.5 ao (1 ao = 
0.529A). The part Vdi sp (r) includes Cg, Cgj C10 terms and 
retardation, while V exc h{r) is given by the analytic form 
1 j r 7/2a—\ exp(— 2ar) derived by Smirnov and Chibisov 
M, in which \o? is the ionization potential of the Rb 
atom in atomic units (au). 

The triplet potential is subject to a larger uncertainty. 
For its short range part an ab initio potential is usu- 
ally taken. To get rid of this model dependence, we use 
the accumulated phase method [ [l0| : the 'history' of the 
atom-atom motion is summarized by a boundary con- 
dition at an interatomic distance r , in the form of the 
phase 4>t(E, of the oscillating triplet radial wave func- 
tion i/j depending on energy E and angular momentum I. 
Specifying (f>r(E,l) is equivalent to giving the logarith- 
mic derivative ip 1 '/ip at r = ro- In all of our previous work 
we neglected the singlet-triplet mixing by the hypcrfinc 
interaction Vh / of the nuclear and electronic spins in the 
range r < ro, in order to deal with pure singlet and triplet 
radial waves until the boundary. Here, however, we in- 
troduce a new variant that allows us to choose a larger ro 
than would otherwise be possible: we include the adia- 
batic mixing by Vh / in the two-atom spin states but still 
neglect its influence on the radial wave functions to avoid 
dependence on the history other than via the pure triplet 
phase. Model calculations show that in this form the 
scattering calculations have the required accuracy for ro 
values up to 16 ao- The experimental data for either ul- 
tracold or weakly bound atoms that we analyze comprise 
a small E and I range near E = I = 0. In this range a first 
order Taylor expansion 4>t{E,1) = + E(j)^, + 1(1 + 1) (j) l T 
is adequate, which reduces the information contained in 
Vr(r) for r < ro to three phase parameters only. In prin- 
ciple, these would be needed for both the 85 Rb and 87 Rb 
systems. However, since we expect Born-Oppcnhcimcr 
breakdown effects to be negligible also for the triplet 
channel in the distance range r < ro, we use mass scaling 
to express <j>^ , <j> l T for 85 Rb in terms of the three phase 
parameters for 87 Rb. Beyond ro we construct Vt(t) from 
V s (r) by adding 2V exch {r). 

Applying this method we carry out a full quantum 
scattering calculation for a set of eight experimentally 
measured quantities. This set consists of five quanti- 
ties for 87 Rb and three for 85 Rb. The 87 Rb data are 
the four bound state energies and the ratio of scatter- 
ing lengths ai_i/d2i = 1.062(12) for atomic scattering 
in condensates of 87 Rb atoms in the hyperfine states 
(f,m f ) = (1,-1) and (2,1) Q. For 85 Rb we include 
the Feshbach resonance fields Bq and B' , as well as the 
energy 0.7(1) mK of the g-wave shape resonance observed 
in the scattering of a pair of cold atoms in the total spin 
S 
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With a least-squares search routine we determine 
optimal values for the parameters Cq, Cg, J, 0! r ( 87 Rb), 
<fi j, ( 87 Rb), <j) l T ( 87 Rb). C10 is kept fixed at the value calcu- 



lated by Marinescu et al. |13|], but the effect of ±10% vari- 
ations around this value and an estimated upper bound 
for the influence of higher dispersion terms are included 
in the final error bars. Column A summarizes the main 
results of the calculations. We find a value for Cq in 
agreement with the theoretical value 4691(23) obtained 
by Derevianko et al. . The C$ value agrees with that 
calculated by Marinescu et al. [jl3j . To our knowledge this 
is the first experimental determination of Cs- Our analy- 
sis also yields the first experimental value of the strength 
of V exc h for cold atoms. The coefficient J agrees with 
the most recent theoretical value in Ref. |L5]]. Table || 
also gives the values of the pure singlet and triplet scat- 
tering lengths for both 85 Rb and 87 Rb, following from 
Cg, Cs, J, 0j.( Rb), as well as the fractional vibrational 
quantum numbers at dissociation vd and the numbers of 
bound states n&. The reduced minimum \ 2 value is 0.5. 

The foregoing makes clear that a major step forward 
has been made possible by the new experiments, two of 
which make use of a Bose-Einstein condensate. This is a 
firm basis for making a variety of interesting predictions. 
As a first example we predict the 87 Rb / = 1 spinor con- 
densate to be ferromagnetic, i.e., it is favorable for two 
/ = 1 atoms to have their spins parallel, because the 
mean field interaction is more repulsive for total F = 
than for F = 2: the calculated scattering lengths are 
a F=2 = +100.4(l)ao and 0^=0 = +101.8(2)a . In a 
recent preprint Klausen et al. independently come 
to this conclusion of a ferromagnetic spinor condensate 
by calculating the scattering lengths for several assumed 
numbers of triplet bound states. 

We are also able to predict collisional frequency shifts 
in an 87 Rb fountain clock for arbitrary choices of partial 
densities of atomic hyperfine states. Table [n] compares 
our calculated fractional frequency shifts normalized to 
total atom density n for two recent experiments |0,E8|. 
We find good agreement with the three measured shifts. 

For various applications there is widespread interest 
for predictions of magnetic field values at which Fesh- 
bach resonances are to be expected in the scattering of 
two 87 Rb atoms in the (/, rrif) — (1, +1) state. With our 
interaction parameters we expect them at the four res- 
onance field values Bo given in Table III together with 
the widths A. The Bq values are to be compared with 
the values 383, 643, 850, and 1018 G predicted in 1997 
M . It is interesting that the broadest resonance at 1004 
G shows a doublet structure [^9) . 

Figure [t] shows Feshbach resonances that we predict 
to occur in the mixed spin channels (2,+l)+(l,-l) and 
(2,-l)+(l,+l) at easily accessible field values of 1.9 and 
9.1G, respectively. The graphs show the predicted field- 
dependent scattering lengths a(B), which are complex 
functions due to the presence of exothermal inelastic de- 
cay channels. The generalized analytic expression for the 
field dependence in this case is |19|: 
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with A e i and Ai ne i the (in)elastic resonance widths and 
4>r a resonance phase constant, arising due to inelas- 
ticity. Note that the real part of the scattering length 
does not go through infinity. It turns out that the 1.9G 
resonance is an I = 2 resonance, which couples via the 
spin-spin interaction V ss to the s-wave incident channel. 
Actually, this resonance is the 'hyperfine analog state' 
p9| in the (2,+l)+(l,-l) scattering channel of the d-wave 
shape resonance of 87 Rb occurring in the spin-stretched 
(2,+2)+(2,+2) spin channel |l2|] , i.e., a state with essen- 
tially the same spatial dependence and differing only in 
its hyperfine spin structure |2Cj . It is located at a compa- 
rable low energy above threshold. In a similar way the I 
= resonance at 9.1G is the hyperfine analog state of two 
of the I = bound states observed || at roughly 25MHz 
below threshold in the (2,+2) + (2,+2) and (1,-1)+(1,-1) 
channels, belonging to the same rotational band as the 
d-wave shape resonance. They might play a role in the 
damping of coherent spin oscillations of the type which 



are being observed in experiments at JILA 21 




2 3 
B (G) 





1 1 1 1 1 1 




f >~ 


< B > ,,,,,, 


, , 1 , , , , 



B (G) 

FIG. 1. (A) Real (solid line) and imaginary (dashed 
line) parts of the scattering length a(B) in s-wave 87 Rb 
(2,+l) + (l,-l) mixed spin scattering channel, showing the 
presence of a Feshbach resonance at 1.9G. Imaginary part 
is proportional to summed rate coefficient G for decay into 
all open channels. (B) Same for Feshbach resonance in 
(2,-l) + (l,+l) channel at 9.1G. 

Until now we assumed the expression for Vdisp to be 
valid for interatomic distances larger than rs- We now 
extend Vdisp with Cn and C12 terms and assume it to 
be valid also between 18 do and r$- This leads us to a 
more ambitious approach that allows us to determine Cio 
and Cn as two more free parameters in the least squares 



search: we take into account the additional constraint 
arising from the equality Vs(r) = V disp (r) - V exch (r) 
by imposing this equality at five r points as additional 
'experimental data' with a standard deviation of 0.5%. 
We thus effectively include the bound states of Ref. Q] 
with outer turning points in the range considered. In 
the search we take C12 equal to the theoretical value 
11.9xl0 9 ofRef. Q. In column B of Table | the resulting 
optimal parameter values are given together with error 
bars based on a 25 % uncertainty in Cn- We find a value 
for Cio differing from the theoretical value 7.665 x 10 7 
au |l3|] by only 1.8%. While the above attractive C„ 
terms with even n arise from the interatomic multipolc- 
multipolc interaction in second order, a Cn term is ex- 
pected [^H as a repulsive third order dispersion term aris- 
ing from the mutual dipole excitation and deexcitation 
of the atoms with an intermediate quadrupole transition 
between excited states in each of the atoms. Note that 
the ratio C11/C12 = —0.072 is comparable to the rigor- 
ous value -0.028 for H atoms and the ab initio ratio 
-0.041 for Cs atoms |^5|. The remaining residue of the fit, 
concentrated at the smallest radii in the radial interval 
may well be due to the summed contributions of further 
(attractive and repulsive) dispersion terms beyond the 
C12 contribution plus correction terms to the Smirnov- 
Chibisov exchange expression. Note that the values of the 
lower dispersion coefficient are dominated by the close- 
to-threshold measurements, whereas the higher ones are 
determined primarily by the Seto potential in the middle 
range ro < r < rs- We expect that experiment will prove 
the value of this more ambitious approach. 

For completeness we point out that a weak contribu- 
tion to the total atom-atom force is still missing in the 
above picture: the interatomic spin-spin interaction V ss . 
One component of V ss is the well-known magnetic dipole 
interaction between the valence electron spins of the in- 
teracting atoms. An additional contribution, which arises 
from the electronic spin-orbit coupling as a second-order 
effect, has been experimentally determined for the first 
time for rubidium atoms by Freeland et al. 0. Calcu- 
lation shows that V ss has a negligible influence on the 
previous analysis. 

In summary, combining the results of three recent high- 
precision experiments we have come close to a complete 
and model independent specification of the interaction 
properties of cold rubidium atoms. Wc have determined 
the van der Waals coefficients Cq, Cs, Cio, Cn, and 
the strength J of the exchange interaction. We have 
thus reached a consistent picture of the interactions, with 
which it is possible to predict essentially all parameters 
needed for a complete description of a rubidium Bose- 
Einstein condensate or thermal gas of any isotope in an 
arbitrary spin state. New experimental data, in partic- 
ular on the Feshbach resonances, will undoubtedly be 
helpful to confirm the above consistent picture and to 
further narrow down the error limits. Wc believe that 
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our approach sets an example for similar experimental 
and theoretical work for other (combinations of) atomic 
species. From a theoretical point of view, it is fascinat- 
ing that it is possible to achieve a level of precision for 
the interaction properties approaching that for hydrogen 
atoms, based on a combination of experimental results 
and a sound framework of collision physics. Additional 
details and their relevance for future experiments will be 
the subject of a future publication p9[ . 
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Quantity 


A 


B 


C e /10 a 


4.703(9) 


4.698(4) 


C 8 /10 5 


5.79(49) 


6.09(7) 


Cio/10 7 


7.665(Ref. Q) 


7.80(6) 


Cn/10 9 




-0.86(17) 


C12/IO 9 
J.10 2 




11.9(Ref. @) 


0.45(6) 


0.42(2) 


a T ( 87 Rb) 


+98.98(4) 


+98.99(2) 


as ( 87 Rb) 


+90.4(2) 


+90.0(2) 


a T ( 85 Rb) 


-388(3) 


-387(1) 


a s ( 85 Rb) 


+2795t£g 


+2400t? 7 ° 


VDT, 71t,T( 87 Rb) 


40.4215(3), 41 


40.4214(2), 41 


VDS,n bS ( 87 Rb) 


124.455(1), 125 


124.456(1), 125 


vdt, n bT ( 85 Rb) 


39.9471(2), 40 


39.9470(1), 40 


vds, nb S ( 85 Rb) 


123.009(1), 124 


123.011(1), 124 



TABLE I. Interaction parameters derived from experi- 
ments without (column A) and including (column B) the re- 
quirement Vs = Vdis P ~ Vexch for ro < r < rs- 



Predicted clock frequency shifts 



(10- 24 cm 3 ) 


Ref. 


Present theory (10 24 cm 3 ) 


-56± 8 t [ 
-50(10)± 22 [ 
-60(16)1 29 | 


17 

H 

18 


-72.5 ± 3.3 
-32.8 ± 0.7 
-41.5 ± 2.9 



TABLE II. Predictions of collisional frequency shifts for 
the 87 Rb fountain clock, compared to two recent experiments. 



Predicted 87 Rb (1,+1) Feshbach resonances 



B (G) 


403(2) 


680(2) 


899(4) 


1004(3) 


A(mG) 


< 1 


15 


< 5 


216 



TABLE III. Resonance fields B and widths A. 
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